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Abstract — A  new  approach  to  MR  angiography,  the  NMR 
diffractive  imaging  technique,  has  been  investigated.  The 
expression  for  NMR  signals  obtained  in  the  NMR  diffrac¬ 
tive  imaging  technique  is  similar  to  the  equation  for  Fres¬ 
nel  diffraction  in  light  waves  or  sound  waves.  Therefore, 
it  is  possible  to  reconstruct  three-dimensional  images  by 
converting  the  NMR  signal  to  hologram  and  using  an  op¬ 
tical  image  processing  system  from  the  data  scanned  two- 
dimensionally.  Moreover,  an  image  focusing  on  optional 
plane  in  the  depth  direction  can  be  reconstructed  numer¬ 
ically  from  the  data  by  changing  an  distance  parameter  in 
the  reconstruction  step.  Experiments  were  performed  using 
an  ultra-low-field  MRI  scanner  to  acquire  two-dimensional 
data  in  the  proposed  technique.  Even  though  blurred  image 
outside  the  focal  plane  is  superimposed  on  the  image  in  the 
focal  plane,  the  three-dimensional  distribution  of  the  object 
can  be  recognized.  This  technique  is  expected  to  be  useful 
in  MR  fast  angiography. 

Keywords — MRI,  angiography,  three-dimensional  imaging, 
holography. 

I.  Introduction 

Lately,  a  brain  scanning  examination,  obtaining  the  im¬ 
age  of  angiography  in  the  brain  to  discover  any  existence 
of  abnormalities,  has  been  utilized  more  frequently.  Imag¬ 
ing  angiography  using  an  MRI  scanner,  however,  requires 
considerable  amount  of  time  to  collect  three-dimensional 
data.  Therefore,  shortening  the  time  for  imaging  is  be¬ 
lieved  to  relieve  the  burden  of  a  study  subject,  followed 
by  improving  the  throughput  because  of  increases  in  the 
number  of  subjects  to  be  checked  per  hour.  With  poten¬ 
tial  possibilities,  shortening  the  time  for  imaging  seems  to 
produce  great  effectiveness  in  the  future. 

In  aiming  to  shorten  the  examination  time,  we  studied 
utilization  of  the  NMR  Fresnel  diffractive  imaging  tech¬ 
nique,  which  enables  us  to  obtain  the  images  containing 
depth  information  by  using  the  collected  two-dimensional 
data.  Although  NMR  phenomena  is  not  a  wave  motion 
itself,  it  can  take  the  similar  form  to  the  description  of 
the  wavefront  in  the  Fresnel  diffractive  region  by  using 
non-linear  field  gradients.  We  have  been  studying  the 
NMR  imaging  technique  in  which  NMR  signal  has  the  form 
of  Fresnel  integral  equation  of  two-dimensional  object  by 
scanning  a  nonlinear  quadratic  filed  gradient  over  the  imag- 
ing  region[l],  [2], 

By  developing  the  Fresnel  transform  imaging  technique 
to  three-dimensional  one,  it  will  be  possible  to  obtain  the 
NMR  images  containing  depth  information  of  the  object  by 
the  signal  scanned  two-dimensionally  as  the  hologram  re¬ 
produces  the  three-dimensional  images  in  the  optical  holog¬ 
raphy. 

Diffractive  tomography  is  a  technique  for  obtaining  sub¬ 
jective  two-dimensional  and  three-dimensional  distribu¬ 


tions  based  on  the  measurement  of  field  waves.  This  tomog¬ 
raphy  has  been  applied  to  diffractive  fields  such  as  optical 
microscopes  [3],  [4]  and  techniques  with  supersonic  waves  [5], 
as  suggested  by  Wolf  concerning  lighting  problems [6]. 

In  our  study,  aiming  principally  to  image  vascular  tracts, 
we  designed  a  coil  system,  which  can  obtain  the  same  signal 
equations  as  those  of  three-dimensional  Fresnel  diffractive 
waves,  and  then  experimented  on  its  use. 

According  to  the  technique  used,  we  can  obtain  the 
three-dimensional  distribution  of  the  subject  in  the  time 
of  two-dimensional  acquisition,  and  we  might  be  able  to 
judge  any  existence  of  abnormalities  rapidly. 

II.  NMR  SIGNAL  IN  THE  FRESNEL  DIFFRACTIVE 
IMAGING  TECHNIQUE 

Holography  is  a  technique  proposed  by  D. Gabor  to 
record  an  object  light  wave  scattered  from  the  object  on 
photosensitive  materials  by  superimposing  a  reference  light 
on  the  object  light  wave.  Although  wavefronts  in  two- 
dimensional  distribution  of  the  object  are  recorded,  illu¬ 
minating  light  wave  on  the  hologram  can  reproduce  the 
scattered  wavefront  from  the  object  which  can  reconstruct 
natural  three-dimensional  images. 

A.  Wavefront  equation  in  the  Fresnel  diffraction  region 

The  diffractive  wavefront  u(xiJyi)  \z=z.  scattered  from 
the  object  g(xo,yo,  zo)  on  the  screen  located  at  z  —  Z{  in 
the  Fresnel  region  is  written  as  Eq.(l)  [8] 

u{xi,Vi) \z=z.  =  ZQj  exp  {jk(zi  -  z0)} 

OO 

I / /g(a?o^o’0°Ux 2(1 -^o)^0^ 

—  OO 

where  A  is  the  wavelength  of  light  source,  k  is  an  wave  num¬ 
ber  and  Zi  represents  the  distance  from  the  center  of  coor¬ 
dinates  (xo^yo^zo)  to  the  screen.  Equation(l)  is  rewritten 
as  follows  by  convolution  equation  : 

u(xuyi)\z=Zi  =  g{xi,yi,Zi)  *  f(xi,yi,Zi),  (2) 

where  f(xi,yi,Zi)  is  written  as  follows: 

f(xi ,  yu  exp  {jk  [z{  +  (x2  +  ?/•  )/2^]  }  .  (3) 

The  input  light  distribution,  g(xo,  yo,  ^o),  has  the  equation 
convoluted  by  quadratic  phase  term  f(xi,yi,Zi)  which  is 
the  function  of  Z{.  The  function  of  f(xi,yi,Zi)  is  called 
“Point  Spread  Function”. 
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B.  NMR  signal  equation 

In  this  section,  we  examine  possibilities  of  an  NMR  sig¬ 
nal  equation  to  obtain  the  similar  equations  as  those  of 
the  three-dimensional  diffracted  wavefront  of  Eq.(l).  Since 
producing  the  same  form  signals  as  those  of  Eq.(l)  makes 
magnetic  field  designing  quite  difficult,  we  planned  to  find 
ways  to  make  magnetic  fields  produced  as  easy  as  possi¬ 
ble  and  approximate  the  function  k/2(zi~  z$)  of  quadratic 
phase  term  in  Eq.(l)  into  a  first  order  equation  as  2  coordi¬ 
nate.  Based  on  the  condition  above,  we  use  a  field  gradient, 
whose  field  intensity  changes  in  a  quadratic  form  on  x  —  y 
plane,  and  whose  coefficient  of  the  quadratic  field  gradient 
varies  in  the  2  direction  by  a,  as  shown  in  Eq.(4) 

A B  =  6(1  +  az)  { (x1  -  x )2  +  (y1  -  y)2 }  ,  (4) 

where  b  is  a  coefficient  of  a  quadratic  field  gradient  at  2  =  0 
and  coordinates  (V,  y')  are  the  center  of  this  quadratic  field 
gradient  which  can  be  set  to  optional  places  by  the  field 
given  from  outside.  Therefore,  Eq.(5)  is  the  target  equation 
of  NMR  signals,  which  can  be  obtained  by  scanning  this 
quadratic  field  gradient  on  xy  plane 


00 


—00 


•  exp  {—  j^br(l+az)  \(xr  —  x)2  +  (y'  —  y)2]  }  dxdydz,  (5) 

where  p(x,y)  represents  the  spin  density  distribution  in 
the  subject,  7  is  the  magnetogyric  ratio,  r  represents  its 
impressed  time,  and  P  is  a  constant. 

In  order  to  obtain  the  signal  like  the  form  of  Eq.(5),  the 
coefficient  of  the  field  must  be  varied  in  the  2  direction 
not  only  the  scanning  field  but  also  the  field  applied  in  the 
time  reading  direction.  Therefore,  magnetic  field  equations 
should  all  be  varied  into  Eqs.(6)  to  (9), 


subject,  a  quadratic  field  gradient  (a)  and  scanning  field  (b) 
are  applied  in  the  phase  encoding  direction  only  for  r 
of  time.  Scanning  field  is  used  to  scan  the  center  of 
a  quadratic  field,  so  when  the  amount  of  scanning  is 
set  to  be  x',  the  quadratic  field  after  scanning  is  A b  = 
6(1  +  az)  {(V  —  x)2  +  y2}.  Next,  the  sweeping  field(c) 
and  the  field  gradient (d)  are  applied  at  the  same  time.  The 
echo  signal  obtained  by  the  method  above  can  be  written 
as  Eq.(10)  by  setting  the  origin  of  the  time  at  the  center  of 
the  echo  signal 


00 

y')  =J J  j  p(x,  y,  z)  exp  {  —j'ybr  (1 +az)  [(x1  -x)2+y2]} 
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When  the  parameters  are  set  to  Eq.(ll),  and  variables  are 
also  transformed,  Eq.(10)  can  be  rewritten  as  the  following 
equation  (12),  and  the  target  equation  of  Eq.(5)  is  obtained 


Gy  boy  f 
=  y 


Gyt 
26r  ’ 


(11) 


OO 

j(x',y')=pJ  J  J  p(x,y,z) 


(a) quadratic  field:  Ab(x,  y,  z)  =  b(l  +  az)  (x2 +y2)  ;  (6) 

(b)  scanning  field:  A b0x  (z)=yj  1  +  azb0x ;  (7) 

(c) sweeping  field:  Aboy(z)  =  (1  +  az)  boy ; ,  (8) 

(d) field  gradient:  A Gyz(y,  z)  =  (1  +  az)  Gyy.  (9) 

Figurel  shows  the  pulse  sequence  for  the  NMR  diffrac¬ 
tive  imaging  technique.  After  excitation  of  spins  in  the 
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Fig.  1.  Pulse  sequence  for  NMR  Fresnel  diffractive  imaging  tech¬ 
nique. 


•  exp  {—jjbr  (1  +  az)  [(xr  —  x)2  +  (y'  —  y)2]  }  dxdydz. (12) 

Where  P  is  set  as  P  =  exp  {—  jjboy/2to2}  and 
jbr(l  +  az),  equivalent  to  a  distance  parameter,  became 
the  same  form  as  written  in  Eq.(5).  There  are  two  methods 
for  reconstructing  images  from  Fresnel  integral  equation. 
One  is  a  technique  using  the  inverse  Fourier  transforma¬ 
tion  once  after  multiplying  the  quadratic  phase  term.  The 
other  method  solve  the  convolution  integral  by  the  inverse 
filteringfl].  In  this  paper,  we  will  discuss  the  method  recon¬ 
structing  images  by  the  inverse  filtering  technique,  in  which 
the  pixel  width  is  not  changed  by  the  imaging  parameters 
depending  on  the  focal  plane. 

Let  Txy  be  the  Fourier  transformation  with  respect  to  x, 
y ,  and  it  is  applied  to  Eq.(12),  then  Eq.(13)  can  be  obtained 
as  follows: 


J  X 
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Fig.  2.  Outlook  of  the  coil  system  for  NMR  Fresnel  diffractive  imag¬ 
ing  technique. 


where  R(kx,ky,  z)  denotes  the  Fourier  transform  of 
p(x,  y ,  z)  with  respect  to  x  and  y  coordinates.  By  arranging 
Eq.(13),  Eq.(14)  is  obtained 


—  oo 

R(t„  k”  Z)  exp  f  47t(l++L-)  } dz-  (14) 

Image  focused  on  z'  plane  is  obtained  by  multiplying  the 
inverse  function  of  modulation  transfer  function  on  z'  plane 
to  Eq.(14)  using  the  distance  parameter  of  rybr(l-\-az'),  and 
taking  the  inverse  Fourier  transform  of  it.  Let  p(x,,y',z,y 
denote  the  image  focused  on  z'  plane  which  can  be  obtained 
by  the  following  equation: 
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From  Eq.(16),  we  can  obtain  the  image  focused  on  optional 
plane  by  giving  the  focal  plane  coordinate  z'  in  the  distance 
parameter  of  jbr^  +  az)  in  the  reconstruction  procedure. 


III.  Experiments 

To  support  this  imaging  technique,  a  coil  system  which 
generate  the  fields  written  by  Eqs.(6)  to  (9)  was  designed 
and  fabricated.  Fig. 2  show  the  outlook  of  the  coil  system. 
Experiments  were  performed  using  an  ultra-low-field  MRI 
scanner  which  generates  the  static  magnetic  field  of  B0  = 
0.0183  T  by  the  solenoid  coil  (resonant  frequency  /0= 779 
kHz).  The  parameters  of  experiments  are  as  follows  :  7 br= 
1.49  rad/cm2,  the  repetition  time  for  the  pulse  sequence 
TR  =300msec,  Axf=Ayf  are  0.2  cm,  and  the  data  matrix 
of  the  NMR  signal  is  set  to  64  x  64. 


Fig.  3.  Reconstructed  images  of  4-pole  phantom  focused  on  different 
depth  by  changing  the  imaging  parameter. 


IV.  Image  Reconstruction 

A.  Adjusting  focal  plane  experiments 

We  conducted  the  experiment  in  order  to  examine 
whether  images  at  optional  focal  planes  can  be  obtained 
from  two-dimensional  scanned  signals  by  parameters  used 
in  the  numerical  reconstruction  procedure.  Fig. 3 (a)  shows 
the  results  of  experiments  using  the  phantom,  having  four 
water  poles  placed  in  the  same  intervals  in  the  x  direction 
and  the  2  direction.  As  the  distance  parameters  used  in  nu¬ 
merical  reconstruction  are  followed  by  jbr  (1  +  az),  com¬ 
puterized  images  was  reconstructed  by  being  varied  from 
1.44  to  1.56  rad/cm2  in  0.44  steps,  as  shown  in  Fig. 3.  Fig. 3 
shows  the  images  focusing  on  a  water  pole  located  in  the 
place,  z  —  Q  mm,  and  (b)  to  (d)  show  the  reconstructed 
images  focusing  on  each  water  pole  after  focal  plane  moves 
by  —0.06cm.  The  results  show  that  the  water  pole  on  focal 
plane  is  imaged  most  clearly.  On  the  other  hand  image 
of  water  pole  become  obscure  as  a  focal  plane  is  moving 
gradually  farther  from  the  plane  where  water  pole  is  lo¬ 
cated.  The  results  of  this  experiment  indicates  that,  even 
in  collecting  two-dimensional  scanned  signal,  NMR  signal 
in  the  NMR  diffractive  imaging  technique  can  reconstruct 
images  focusing  on  optional  plane  according  to  the  choice 
of  distance  parameters  depending  on  2  coordinate. 

B.  Imaging  experiments  with  3D  tube 

We  performed  the  experiment  using  a  vascular  tract 
model.  Phantom  shown  in  Fig. 4  was  used  for  this  experi¬ 
ment  after  being  devised  with  a  4  mm  diameter  tube  includ¬ 
ing  water.  Numerically  reconstructed  images  obtained  by 
adjusting  the  focal  plane  z'  are  shown  in  Fig. 4(a)- (d)  whose 
focal  plane  is  shifted  every  —8mm  from  a  point  z  =  12 
mm.  Fig. 4(a)  and  (b)  show  the  images  each  focusing  on 
the  tubes  on  the  right  and  left  sides  of  a  small-sized  loop. 
Fig. 4(c)  and  (d)  show  the  images  each  focusing  on  the  tubes 
on  the  right  and  left  sides  of  a  large-sized  loop.  Although 
these  images  include  unnecessary  image  components  out  of 
focus,  the  three-dimensional  distributions  of  phantom  are 
understood  by  three-dimensional  information. 

Using  a  proposed  diffractive  imaging  technique,  there  are 
some  cases  where  images  are  emphasized  by  the  superimpo¬ 
sition  of  blurred  image  components  outside  the  focal  plane 
on  the  focal  plane  one.  The  area  indicated  as  A  in  Fig. 4(d) 
is  the  example  of  that  obscure  image.  In  this  case,  it  may 
not  be  possible  to  distinguish  normal  parts  from  diseased 
parts.  In  order  to  reduce  this  interference  problem,  we  had 
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Fig.  5.  Front  view  images  of  3D  tube  phantom  focused  on  different 
depth  by  changing  the  imaging  parameter. 


a  hologram  which  can  produce  image  in  the  coherent  op¬ 
tical  system  is  produced  easily  by  transforming  the  signal 
distribution  into  the  gray-scale  pattern  [9].  In  principle,  it 
is  possible  to  reconstruct  natural  three-dimensional  images 
using  the  holographic  technique.  We  will  investigate  the 
holographic  reconstruction  in  the  next  step. 


Fig.  4.  Diagonal  view  images  of  3D  tube  phantom  focused  on  different 
depth  by  changing  the  imaging  parameter. 

performed  the  experiments  for  imaging  from  different  an¬ 
gles  to  change  the  way  of  interference.  Fig. 5 (a)- (d)  show 
the  results  of  an  imaging  experience  by  revolved  phantom. 
Since  the  way  of  interference  is  changed,  it  ease  the  judge¬ 
ment  whether  parts  with  high  intensity  is  generated  by 
interference  or  by  other  reasons. 

V.  Discussion 

In  the  experiment,  we  could  obtain  images  focusing  on  an 
optional  depth  by  the  signals  scanned  two-dimensionally. 
We  can  recognize  that  an  object  is  located  around  focal 
plane  by  the  fact  that  the  image  is  clearly  reconstructed. 
However,  there  are  some  cases  where  we  have  found  it  dif¬ 
ficult  to  judge  this  clearly.  In  this  case,  we  can  recognize 
three-dimensional  distributions  of  an  object  more  easily, 
when  we  can  observe  images  as  animation  by  continuously 
displaying  the  images  focused  on  each  2  coordinate.  This 
is  considered  to  be  because  we  can  spatially  recognize  the 
obscure  amount  of  images  located  in  the  front  and  rear 
parts  of  focuses  by  observing  continuous  images  in  the  2 
direction. 

However,  given  the  considerable  shortage  of  information 
about  the  subject  on  the  use  of  this  technique,  we  will 
have  to  confront  problems  of  judging  whether  an  observed 
place  is  abnormal  because  the  image  intensity  on  the  plane 
will  become  greater  when  vascular  tracts  are  crossing  in 
the  depth  direction.  In  this  case,  although  it  doubles  the 
amount  of  imaging  time,  we  might  solve  the  problem  by 
observing  the  subject  from  different  angles.  There  are  two 
ways  to  move  a  visual  point:  one  is  to  move  an  object  and 
the  other  is  to  create  an  equivalent  effect  as  moving  a  visual 
point  by  adding  magnetic  fields  in  the  imaging  experiment. 
We  consider  that  the  latter  technique  -  moving  visual  points 
by  fields  -  is  the  more  practical. 

Since  the  description  of  the  NMR  signal  has  the  similar 
form  to  the  equation  of  the  wavefront  in  the  Fresnel  region, 


VI.  Conclusion 

A  new  approach  to  MR  angiography,  the  NMR  Fresnel 
diffractive  imaging  technique,  which  can  obtain  the  signal 
similar  to  the  equation  of  the  Fresnel  diffractive  images,  is 
proposed.  To  support  this  imaging  technique,  a  coil  sys¬ 
tem  composed  of  six  coils  was  designed  and  imaging  exper¬ 
iments  using  the  water  tube  model  were  performed.  The 
results  shows  that  images  focusing  on  optional  plane  in  the 
depth  direction  can  be  reconstructed  from  data  scanned 
two-dimensionally.  Even  though  blurred  images  outside 
the  focal  plane  are  superimposed  on  the  image  in  the  focal 
plane,  the  three-dimensional  distribution  of  the  object  can 
be  recognized  by  moving  the  focal  plane  in  the  depth  di¬ 
rection.  To  attain  supplemental  information  for  the  object, 
acquiring  images  from  different  angles  is  helpful  for  recog¬ 
nizing  the  spatial  distribution  of  the  object  more  precisely. 
The  proposed  imaging  technique  is  expected  to  be  useful 
in  MR  fast  angiography. 
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